Scaling of Entropic Shear Rigidity 
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The scaling of shear modulus near the gelation/vulcanization transition is explored heuristically 
and analytically. It is found that in a dense melt the effective chains of the infinite cluster have sizes 
that scale sub-linearly with their contour length. Consequently, each chain contributes k^T to the 
rigidity, which leads to a shear modulus exponent dv. In contrast, in phantom elastic networks the 
scaling is linear in the contour length, yielding an exponent identical to that of the random resistor 
network conductivity, as predicted by de Gennes. For non-dense systems, the exponent should cross 
over to dv when the percolation correlation length is much larger than the density-fluctuation length. 
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Introduction — Gelation and vulcanization are continu- 
ous phase transitions from liquids to random solids, 
caused by the introduction of chemical (i.e. permanent) 
crosslinks; see, e.g., Ref. [l|. It is by now well established 
that the geometrical aspects of these transitions are cor- 
rectly described by percolation theory 0, 0, 13 ■ Qualita- 
tively speaking, both transitions — gelation/vulcanization 
and percolation — are controlled by the emergence and 
structure of an infinite cluster at the critical point. 

The elastic and thermodynamic properties of gels near 
the critical point are not as well established. Of these, 
the most controversial is the scaling behavior of the 
static shear modulus /i, which is defined in the follow- 
ing way. Consider a spatially homogeneous, volume- 
preserving shear deformation A [J], with det A = 1. Un- 
der such a deformation, the increase in the free energy of 
a gel is, to leading order, given by Q 




5F = V^l{TY^-d). 



g = A -A, 



(1) 



where g is the metric tensor and d is the spatial dimen- 
sionality. The shear modulus is expected, on general 
grounds, to obey a scaling law near the critical point. 



M = Mo0( 



(2) 



where the control parameter r measures (minus the) de- 
viation of the cross- link density from criticality (i.e. r < 
in the solid phase), and Q is the unit step function. 

Values reported for the exponent /, either from ex- 
periments or computer simulations, are rather scattered, 
and seem to suggest four different universality classes [GJ . 
For systems with entropic elasticity, which is the focus 
of this Letter, values of / usually fall into one of two 
classes. Firstly, most numerical simulations 7] involv- 
ing phantom networks, as well as many gelation exper- 
iments, suggest that / has the same value as the con- 
ductivity exponent i (w 1.9 in three dimensions) of ran- 
dom resistor networks (RRN), supporting a conjecture 
of de Gennes' [jj]. A second class of experiments, as well 
as some simulations, support the scaling result / = dv 



FIG. 1: The infinite cluster is a network of effective chains: 
(A) In a dense melt each effective chain is one single thermal 
blob with typical size ^perc, and contributes k^T to the shear 
rigidity. (B) In a phantom network each chain is sequence 
of "tension blobs" (dashed circles), each of size ^ton. Each 
tension blob contributes fceT to the shear rigidity. 



(« 2.6 for three dimensions), where v is the percolation 
correlation-length exponent, as proposed in Ref. jg. The 
argument for the latter exponent is in the spirit of the 
classical theory of rubber elasticity (see, e.g., Ref. pi), 
in which the elastic modulus acquires a contribution of 
fceT per effective chain. The purpose of this Letter to use 
heuristic and analytical methods to outline a resolution 
of the long-standing apparent contradiction between the 
two plausible arguments, mentioned above, as well as the 
inconsistency across experimental and simulational data. 
Heuristic reasoning — Let us first consider gelation in a 
dense system with strong inter-particle repulsions, and 
let the system be near the critical point. Then the 
characteristic length-scale for density fluctuations ^don 
is much smaller than the percolation correlation length 
Cpcrc (beyond which the infinite cluster is effectively ho- 
mogeneous l_^Now let us invoke the "nodes-links-blobs" 
picture [9|, [lOJ , in which the incipient infinite cluster is 
a network of effective chains, connected to one another 
at effective vertices, called nodes. For our purposes, it 
is adequate to treat the effective chains as quasi-one- 
dimensional objects having some average thickness. The 
end-to-end displacement of these effective chains has a 
certain distribution with a characteristic length-scale, 
which is presumably identical to the percolation corre- 
lation length ^pcrc- In the absence of external stress , w e 
expect these chains to exhibit a type of random walk [ll| , 
owing to thermal fluctuations. Therefore, the end-to-end 



displacements should scale sub-linearly with the contour 
length. In the language of polymer physics, every effec- 
tive chain constitutes a single "thermal blob"; see Fig.Q] 
and Ref. [13j- Under a small shear deformation, each 
chain will be slightly stretched or compressed, thus con- 
tributing fceT to the shear modulus. As there are roughly 
effective chains per unit volume, the shear modulus 



^perc 



should scale as kBT^~^j.^ ~ \r\ , near the critical point. 
We note that the shear modulus has dimensions of en- 
ergy per unit volume; hence, this scaling is also mandated 
by dimensional analysis, provided that ^porc is the only 
important length-scale near the transition. 

Let us now consider typical numerical simulations of 
phantom networks, i.e., systems without inter-particle 
repulsions. The crucial observation, made above, that 
the equilibrium conformations of the effective chains are 
un-stretched random walks (i.e. "thermal blobs") with 
typical size ^porc now breaks down, for the following rea- 
son. In such simulations, the sol fraction (i.e. the finite 
clusters) are usually removed completely, as they do not 
interact with the gel fraction. The resulting gel is very 
sparse, and would tend to collapse so as to maximize 
the entropy. To prevent this, the system size is usu- 
ally fixed during shear deformation. Consequently, long 
effective chains are strongly stretched. The conforma- 
tion of the infinite cluster is such that the net entropic 
force at each node vanishes, i.e., the cluster is statisti- 
cally in mechanical equilibrium. Therefore, the mean 
tension S carried by each chain has the same order of 
magnitude. As shown Fig. ^, 5" defines a length-scale, 
^tcn = k^T / S , beyond which chain conformations are 
dominated by tension, and are thus effectively straight. 
By contrast, within ^ten thermal fiuctuations dominate, 
so that chain conformations are random walks. £tsp is, 
by definition, the typical size of a "tension blob" [IjI for 
a polymer under tension S. Near the critical point, ^tcn 
is much smaller than ^perc, which diverges as Irl^"^. It 
follows that the number of tension blobs on each chain, 
as well as the chain's end-to-end displacement, scales lin- 
early with its contour length in a phantom network. Un- 
der a shear deformation, every tension blob contributes 
k^T to the overall shear rigidity [l3|. Therefore, a typ- 
ical chain, comprising many tension blobs, contributes 
a term to the total shear rigidity that is proportional 
to its contour length. This should be contrasted with 
the case of melt, in which every chain contributes ksT , 
independent of its contour length. In a coarse-grained 
description, we may replace every tension blob by a me- 
chanical spring of natural length zero and force constant 
/cB^'/Cten' without changing the elasticity. The result- 
ing model is a randomly diluted network of mechanical 
springs of zero natural length at zero temperature, which 
can be mapped into the randomly diluted resistance net- 
work model [lj|. In this mapping, the coordinates of 
nodes are mapped into voltages, and the shear modulus 
into the conductivity. Therefore the shear modulus of 



a randomly diluted entropic phantom network is equiva- 
lent to the conductivity of a random resistor network, as 
de Gennes conjectured and many numerical simulations 
have supported. 

Analytical reasoning — We have extended the Landau the- 
ory for the elasticity of vulcanized matter[2| to the case 
of tunable repulsive interactions |l5| . The relevant order 
parameter is the (1 -I- n)-fold rephcated particle density 
distribution, less its expectation value in the liquid phase: 

n{x) = rj(x",x\ . . . ,x") = -Ary-(i+") (3) 

+ E,1i ['^(^° " ^^■)^'5(^' - ^^■)) • • • ('5(^" - ^^■))] 

Here, x is short-hand for (1 -I- n) d-dimensional vectors 
(x°,x^, . . . ,x"); Cj labels the position of the j^^ parti- 
cle in the system, (•) denotes a thermal average over the 
measurement ensemble 16], and [•] denotes the average 
over the cross-linking realization (i.e. the quenched disor- 
der) as well as the thermal fluctuations of the preparation 
ensemble. J7 can be interpreted as giving the joint prob- 
ability density function (p.d.f.) that a particle is located 
at x" in the preparation state and is later found at po- 
sitions x^ , . . . , x" in n independent measurements in the 
measurement state. Q,{x) encodes a great deal of informa- 
tion about random solid state. In particular, the particle 
density fluctuations of the preparation and measurement 
states are respectively given by 

f)"(x)=E5(x~c,)]-7V/V^, (4a) 

f]"(x)=[(^5(x-c,))]-7V/U, (4b) 

n-{^-)^j\['^^^^^^^d^,n{x). (4c) 

In the liquid phase, all particles are delocalized, i.e., 
((5(x" — Cj)) = V^^ for all j. This ensures that the order 
parameter 51(a;) vanishes identically. In the gel phase, 
however, a certain fraction of particles belong to an infi- 
nite cluster and are localized. If such a particle were at 
position x" in the preparation state, it would fluctuate 
around the same point in the measurement state [mod- 
ulo a global translation and rotation that are common to 
the system as a whole). This is captured by a nonzero 
expectation value of the joint p.d.f. ||2J) in the gel phase. 
The Landau free energy functional comprises one part 
concerning localization, and another describing density 
fluctuations. The localization part accounts for the ef- 
fects of cross-links, and is given by 

H^=j dxkn^+':f{y'nf+'^J2{y-nf^^,nA, 

(5) 
where V*^ and V" are, respectively, derivatives with re- 
spect to yp and x" . The parameter r controls the cross- 
link density and drives the transition to the random 
solid state, whereas Kq and K respectively measure the 



stretchability of polymer chains in the preparation and 
measurement states (the larger the iC's, the softer the 
chains). The free energy cost for density fluctuations de- 
pends only on one-replica quantities, f2"(x"), and is 

Hd = ^ j d^^n\^°f + ^Y. j d^-n-{^'^)\ (6) 

a— 1 

where the non-negative parameters Bq and B are, respec- 
tively, the compressibility of the system in the prepara- 
tion and measurement states. A large value of Bq would 
ensure that the system is cross-linked in a state with al- 
most uniform density profile, i.e., vanishingly small fluc- 
tuations in ri". On the other hand, B characterizes the 
repulsive interactions between particles in the measure- 
ment state. In a typical vulcanization experiment on a 
concentrated solution or melt, both B^ and B are large, 
so that the density remains essentially uniform across the 
transition, ft is important, however, to realize that Bq 
and B are separately adjustable (as are Kq and K), e.g., 
via tuning the solvent quality before and after crosslink- 
ing. For gelation, the network formation process com- 
monly lasts for an extended period. As a result, B may 
differ from Bq, even by a large factor (e.g. due to correla- 
tions built up during the course of the chemical reaction) , 
even if all external physical conditions remain unchanged. 

fti strong contrast, in typical numerical simulations of 
phantom systems, all polymers (or particles) are ascribed 
to lattice sites, and then cross-links are randomly intro- 
duced, connecting some neighboring particles. After re- 
moving the sol part, the system is allowed to relax at 
nonzero temperature, with intra-cluster repulsion com- 
pletely ignored. This corresponds to a large positive 
value for Bq but a vanishing value oi B. As we shall 
soon see, it is this qualitative difference between Bq and 
B that is responsible for conductivity-like scaling of the 
shear modulus in phantom networks. 

To study the elastic properties, we consider deforming 
the system after cross-links have been introduced. As 
shown in Fig. [SJ this amounts to making an affine change 
A of the boundaries for measurement replicas (1 through 
n), leaving the preparation replica intact, ft is convenient 
to then make the linear coordinate transformation 

c"->A-c", x"->A-x", forl<a<n, (7) 

which restores the original boundary conditions for the 
n measurement replicas. Under this transformation, Hd 
of Eq. ||nj), is unchanged, provided we redefine B appro- 
priately. As for Hx, Eq. Q, wc find that all terms are 
invariant under the transformation (in the n — > limit), 
except for the term with coefficient K, which becomes 

n d n 

- ^ {v-^f -^ Y ^ ■9'^^' ^ ^"^ ^^" ' (^) 

a — l a, 6—1 a — 1 

where g is the metric tensor defined in Eq. |Q. 




a=l 



a=« 



FIG. 2: Shear deformation (full vs. dotted lines) is applied 
to measurement but not preparation replicas. After the coor- 
dinate transformation, the original boundaries are recovered 
at the cost of introducing a nontrivial metric tensor. 



The effects of fluctuations, both thermal and quenched, 
are studied by averaging over all order-parameter config- 
urations. Then the disordered-averaged physical free en- 
ergy [F] is related to the (l-l-n)-fold rephcated partition 
function ^i+„ by 



[F] 



lim dn In Zi_ 

n-»0 



lima„ln/ Dne~"^+"^, 

n^O J 



where -ffi+n (= Hj-, + Hx) is the full effective[lg Hamil- 
tonian; see Ref. \W\. 

A renormalization-group (RG) analysis of the full 
model, described by Hi+n, will be presented in a fu- 
ture publication |l5|. Here, we shall only present re- 
sults for two hmiting cases: (a) Bq = B = +cx3, i.e., 
vulcanization in an incompressible polymer melt; and 
(b) Bq = -l-oo but B = 0, i.e., phantom networks. We ap- 
ply a momentum-shell RG transformation, thus integrat- 
ing out short length-scale fluctuations recursively, whilst 
rescaling the order-parameter fleld H, and spatial coor- 
dinates such that Kg and K remain unity. The renor- 
malizations of r and v turn out to be independent of B 
and g, i.e., the parameters describing the measurement 
ensemble. To one-loop order we find: 



dr 
U 



6 



2 ' 



dv 



7 



-ev 



(9) 



where e = 6 — d. For d < 6 there is a nontrivial fixed point 
at (r*,u*^) = (e/14, 2e/7). Correspondingly, the critical 
exponents (77, v) are given by (— ^j-e, 5 + ^^)' which agree 
with results from the the e expansion for the percolation 
transition '20|. 

The flow of g determines the elastic properties. Ex- 
pressing g as (?g, where g has unit determinant [so that 
g"^ = det(g)] we find that, regardless of the value of B, 
g does not flow. Therefore, renormalization of the met- 
ric g is completely controlled by its determinant. For 
B = Bq = +00, we flnd the flow equation 



dg/dl^{2/3)v^l~g')g. 



(10) 



For d < 6, v^ -^ v*^ = 2e/7 > 0, hence g flows to 
unity. Qualitatively, this implies a symmetry between the 
preparation and and measurement ensembles. It also sug- 
gests that the correlation length ^porc is the only relevant 
length-scale, in agreement with the preceding heuristic 
argument, that, in a dense melt, each effective chain 



constitutes a single thermal blob of typical size ^porc- 
Near the fixed point, the singular part of the free energy 
(i.e. the elastic free energy) has the scaling form 



/s = |r|'^''^i(g*) = |r|'^''^i(g). 



(11) 



As the shear modulus is given by an appropriate deriva- 
tive of /s with respect to g [cf. Eq. Q], we immediately 
see that it scales as |r|''''. 

For the second case, viz., {Bq,B) — (+00, 0) we find 

dg/dl = {2/3)v'^g — >4:eg/21 (12) 

Now g, and also the metric tensor g, are relevant near 
the percolation fixed point, with a positive crossover ex- 
ponent (f>g of 4e/21, echoing our heuristic argument that 
effective chains are strongly stretched in a phantom net- 
work. In general, the singular part of free energy should 
then have the scaling form 



fs^lrf^Ms/M^')- 



(13) 



For a pure shear with det(g) = 1, this fs must agree with 
Eq. 1^, up to a constant independent of g. Therefore 
the shear-modulus exponent is given hy dv ~ (f>g = 3 — 
^e, which, to the same order in e, is identical to the 
conductivity exponent of a random resistor network |2fll | . 

This equivalence between critical exponents of phan- 
tom elastic networks and random resistor networks 
should in fact hold to all orders in e. To see this ana- 
lytically, we note that setting {Bq,B) = (00, 0) in Hd is 
equivalent to setting Bq = B = 0, together with the hard 
constraint fl'^{x.) = 0, which explicitly excludes configu- 
rations having nonzero density fluctuations in the prepa- 
ration ensemble. The resulting model then becomes for- 
mally identical to the Harris-Lubensky formulation of the 
random resistor network problem [2(J, with the nd (^ 0) 
coordinates associated with the measurement ensembles 
(x^, . . . ,x") mapped onto the D (— > 0) replicated volt- 
ages d, provided the stated limits are taken. Thus, the 
pair of systems are governed by identical RG equations 
and, hence, critical exponents. 

Having established the two limiting cases, B — 00 and 
-B = 0, it is natural to ask which one is the more stable. 
As B, like r, has naive dimension 2, it is always relevant 
near 6 dimensions. Therefore, if we keep Bq large and 
tune B to be small, the shear modulus exponent should 
cross over — from the conductivity one, t, to the incom- 
pressible system one, dv — when \r\ becomes smaller than 
B. This can be readily realized in a numerical simulation, 
if one were to retain the sol part of the system and turn 
on a small repulsion. In principle, this cross-over might 
also be observed in gelation experiments on non-dense 
solutions, provided |r| is sufficiently small. 
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